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We report the use of broadband heterodyne spectroscopy 
to perform continuous measurement of the interaction energy 
-Eint between one atom and a high-finesse optical cavity, dur- 
ing individual transit events of ~ 250 duration. We achieve 
a fractional sensitivity ~ 4x 10^^/\/Hz to variations in E^y^t/h 
within a measurement bandwidth that covers 2.5 decades of 
frequency (1-300 kHz). Our basic procedure is to drop cold 
Cesium atoms into the cavity from a magnetooptic trap while 
monitoring the cavity's complex optical susceptibility with a 
weak probe laser. The instantaneous value of the atom-cavity 
interaction energy, which in turn determines the coupled sys- 
tem's optical susceptibility, depends on both the atomic posi- 
tion and (Zeeman) internal state. Measurements over a wide 
range of atom-cavity detunings reveal the transition from res- 
onant to dispersive coupling, via the transfer of atom-induced 
signals from the amplitude to the phase of light transmit- 
ted through the cavity. By suppressing all sources of excess 
technical noise, we approach a measurement regime in which 
the broadband photocurrent may be interpreted as a classi- 
cal record of conditional quantum evolution in the sense of 
recently-developed quantum trajectory theories. 



Optical cavity quantum electrodynamics (QED) in the 
strong coupling regime Q provides a unique experimental 
paradigm for real-time observation of quantum dynami- 
cal processes at the single-atom level. While spectacular 
advances have certainly been made in the preparation 
and tomography of quantum states of motion for a single 
trapped ion [0J3] , all such experiments have involved the 
accumulation of ensemble-averaged data over many suc- 
cessive realizations of the process being studied. Recent 
studies of single-molecule dynamics have likewise demon- 
strated the "immediate" detection of photochemical Q or 
conformational |^ events, but such experiments presently 
lack the potential that cavity QED provides for observing 
quantum processes on a timescale that makes coherent 
control/intervention a tangible possibility. We wish here 
to look beyond the mere detection of quantum jumps, 
and to focus on the development of a broadband, single- 
shot measurement technique that achieves signal-to-noise 
ratio > 1 over a bandwidth that includes all characterstic 
frequencies of a quantum dynamical process. 

Real-time observation of quantum dynamics in many- 
atom systems has recently become an important theme 
in atomic physics, with notable demonstrations involving 
vibrational excitations of a trapped Bose-Einstein con- 
densate |6| and the decay of coherent oscillations of an 
ensemble of atoms in an optical lattice pUg] . In contrast 



to programs like these, for which the scientific empha- 
sis lies on noninvasive observation of a system's intrinsic 
dynamical processes, experiments in single-atom cavity 
QED hold great potential for enabling precise investiga- 
tions of how measurement backaction alters the dynam- 
ical behavior of a continuously-observed open quantum 
system 

A sophisticated theoretical basis for understanding 
such issues is presently maturing in the form of quan- 
tum trajectory theories [p^|-p^, but significant techni- 
cal challenges remain to be solved before definitive ex- 
periments can be performed in the lab. Our purpose 
in the present work is to report substantial progress to- 
wards surmounting such obstacles in the context of cavity 
QED, and hence towards achieving the essential experi- 
mental capabilities required to perform quantitative tests 
of measurement-based Stochastic Master Equations. We 
ultimately hope to be able to implement some recently- 
proposed "applications" of the continuous observation of 
dissipative quantum dynamics, in fields such as quantum 
measurement 0,|O|, quantum chaos [[8|jl|], and quan- 
tum feedback control ||l|,|^,||j2|. 

This article focuses on a detailed description of our re- 
cent experiments that record the complete time-evolution 
of interaction energy between one atom and a high-finesse 
optical cavity, during individual transit events of ~ 250 
/is duration. With characteristic atom-cavity interac- 
tion energies E-mt/Ti ~ 10 MHz, we achieve measure- 



ment sensitivities Sg ~ 4.5 kHz/vHz over a bandwidth 
that covers the dominant rates of variation in Ei^t (1- 
300 kHz). Unlike typical pump-probe measurements of 
scattering dynamics in real {e.g. diatomic) molecular sys- 
tems p5| , p6[ |, our experiments on the Jaynes-Cummings 
"molecule" yield a continuous time-domain record of the 
atom-cavity coupling during each individual "scattering" 
event (transit). The data clearly illustrate variations 
caused by atomic motion through the spatial structure 
of the cavity eigenmode and/or optical pumping among 
the atomic internal (Zeeman) states. In certain param- 
eter regimes of the detuning and probe power, distinc- 
tive indications of the quantum-mechanical nature of the 
atom-cavity coupling can be seen in the photocurrent 
recorded from just a single atomic transit. For large 
(> 50 MHz) atom-cavity detunings we obtain phase- 
contrast signals induced by individual atomic transits, 
corresponding to a regime of strong but dispersive cou- 
pling. Phase-quadrature measurements of atomic motion 
have been widely discussed in the quantum optics litera- 
ture |^7 H29| , |lT| , p^ , p^ , but the present work provides the 
first experimental demonstration at the single-atom level. 
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Because of the standing-wave spatial structure of the 
cavity eigenmode, and the corresponding rapid varation 
of the atom-cavity coupling strength over sub- wavelength 
distances, our data should in principle display a sensitiv- 
ity of 1.5 X 10"^" m/VHz to atomic displacements along 
the cavity axis. Unfortunately we cannot realize this fig- 
ure as a precision for monitoring the atomic position, as 
we do not presently have any means of separating sig- 
nal variations due to motion through the standing wave 
from "background" contributions due to transverse mo- 
tion or optical pumping. In our concluding section, we 
shall briefly discuss our motivations for further work to 
disambiguate the nature of rapid variations in our data. 



I. BASIC THEORETICAL DESCRIPTION 

In simple terms, our experimental procedure is to drop 
a cloud of cold Cesium atoms from a magnetooptic trap 
(MOT) into a high-finesse optical cavity, while continu- 
ously monitoring the cavity's complex susceptibility with 
a weak probe laser |3^j3l[]. By limiting the number of 
atoms in the initial cloud, we can easily reach an oper- 
ating regime in which atoms transit the cavity only one 
at a time. Using broadband heterodyne detection and a 
high-speed digitizer, we continuously record both the am- 
plitude and phase of the transmitted probe beam during 
^ 50 ms time windows. Each window typically contains 
from zero to five atom transit signals. 

The elementary theoretical description of such a mea- 
surement employs steady-state solutions of the nonse- 
lective Master Equation for a stationary two-level atom 
coupled to a single electromagnetic mode via the Jaynes- 
Cummings interaction Hamiltonian. Atomic center-of- 
mass motion and optical pumping among Zeeman states 
can only be included in this treatment by allowing for 
a time-dependent atom-cavity coupling strength. Al- 
though this type of approach cannot make predictions 
about dynamical variations in the coupling strength, it 
does provides a quantitative basis for interpreting some 
time-independent features of our data. Our use of this 
"adiabatic" model may be justified to a certain extent 
by the separation of timescales that we achieve in optical 
cavity QED with laser-cooled atoms. In optimal cases, 
the atom-cavity coupling strength should vary by as little 
as a factor of lO""* over the system damping time ^ 30 
ns. 

If we let p denote the density operator for the joint 
state of the atom and cavity, the nonselective Master 
equation (in the electric dipole and rotating- wave approx- 
imations) reads: 



Ho,P + 7_L (2(7/9(7^ — (7^(7/3 — pO-^lj) 

+ i^a + + Kc) {2apa) — a) ap — a) ap) , (1) 

= TiAa'I'a + hQa'^a + \/2k^£ (a + a)) + Hint, (2) 

Hint = ngoe^(y'+'')^'"' cos (yfcix) [aa^ + a^a] . (3) 



Here 7j_ is the atomic dipole decay rate, Ka is the cavity 
field decay rate through the input mirror (through which 
the probe laser is injected), Kb is the cavity field decay 
rate through the output mirror, Kc is the cavity field de- 
cay rate due to intracavity scattering/absorption losses, 
A = t'a — t'p is the atom-probe detuning, O = Vc ^ i^p 
is the cavity-probe detuning, and the coupling strength 
go is equal to half the maximum single-photon Rabi fre- 
quency. We treat the atomic center-of-mass coordinates 
X, y, z as c-number parameters, with the x-axis conciding 
with the cavity axis and z parallel to gravity. The Gaus- 
sian waist of our cavity mode is w ~ 45 pm. Note that we 
have written the Master Equation in a frame rotating at 
the drive frequency, so £ is a constant term proportional 
to the complex amplitude of the driving field. 

To find the steady-state density operator as a function 
of driving strength and various detunings, we simply set 
Pss = and solve for pss using linear algebra. The ex- 
pected amplitude and phase of the heterodyne photocur- 
rent may then be computed as p2p^ 

{^het (t)) = vfL^V^iTbexp [i {^Lt + (t>L)] Tr [p,,a] , (4) 

where 77 represents the overall photodetection efficiency 
(including propagation losses between the cavity and 
photodetectors, heterodyne efficiency, and detector quan- 
tum efficiency), /l and 4>l represent the photon flux and 
phase of the (optical) local oscillator, and fi^ is the fre- 
quency of the optical local oscillator relative to the rotat- 
ing frame (cavity driving field). In the experiment, we 
mix iiict (t) with an rf local oscillator at the frequency 
Hl (which ranges between 40 and 190 MHz) and sepa- 
rately record the in-phase and quadrature components of 
the slowly-varying envelope (with an analog bandwidth 
of 300 kHz). 

Note that the presence of an intracavity atom can in- 
fluence the heterodyne photocurrent only via the inter- 
action Hamiltonian _ffi,it oc hg{r), where 

<y'+-')/^\os{kLx) . 
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In a two-level approximation for the atomic internal dy- 
namics, and for a classical treatment of the atomic center- 
of-mass degrees of freedom, all steady-state properties of 
the atom-cavity system are strictly determined by the 
value of g(r) once the parameters {Ka,Kb,j±, A,Q,£) have 
been specifled. This includes the quantity (Hint), which 
represents one possible measure of the "interaction en- 
ergy" between atom and cavity. Given that ^(r) ap- 
pears to be a more fundamental measure of the interac- 
tion strength, however, we have adopted the convention 
Eint = %(r)- 

The atom-cavity interaction can be treated semiclassi- 
cally using the optical bistability state equation (OBSE) 
psf . The OBSE is traditionally written in terms of the 
scaled field variables x and y (not to be confused with 
the atomic coordinates), with the correspondence 
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where the saturation photon number mo is given by . 
For a given driving strength £, the expected intracavity 
field amphtude can be found by inverting the equation 

where we again work in a rotating frame at the drive fre- 
quency, S = {uja ~ ^p) /l± represents the scaled atom- 
probe detuning, cj) = (wc — (jJp) / (ko + Kb + Kc) repre- 
sents the scaled cavity-probe detuning, and the "coop- 
erativity" C is defined by 



2 (Ha + Kb + Kc)7± ' 

Note that the dependence of g (and thereby C and mo) 
on the atomic position and internal state is implicit. A 
semiclassical prediction for the heterodyne photocurrent 
is obtained by substituting Xy/mo for Tr [ps^a] in equa- 
tion (|). 

II. EXPERIMENTAL APPARATUS AND 
PROCEDURES 

Figure |l| provides a general overview of the apparatus, 
indicating the schematic arrangement of various compo- 
nents to be described below. The diode laser setup for 
forming the Cs MOT is not shown. 



A. High-finesse microcavity 

We use a Fabry-Perot high-finesse microcavity 
("physics cavity") consisting of two spherical mirrors 
with 1 m radius of curvature [^4| . The cavity was con- 
structed with a mean length I ~ 107.5 /im, which we 
inferred from the cavity's measured free spectral range 
of 1.395 X 10^^ Hz. The measured I and specified radii 
of curvature geometrically determine the cavity's electro- 
magnetic mode volume for TEMqo modes near 852 nm 
p5| . Together with the dipole decay rate 7_L/27r ~ 2.6 
MHz for the Cs 6P3/2 level this determines our op- 
timal coupling constant go/2TT to be ~ 11 MHz for a± 
transistions (specifically the 651/2 {F = A,mF = ±4) 
6^3/2 (-P" = 5, mj? = ±5)) and ~ 6 MHz for tt transi- 
tions (651/2 (F = 4, mF = 0) ^ 6A/2 {F = b^mp = 0)) 
within the D2 Zeeman manifold ||l|,p7t. 

In order to allow cold atoms to fall into such a short 
cavity, we found it necessary to have the mirror manufac- 
turer reduce the substrate diameters from the standard 
value of 7.75 mm down to 3 mm |Q. This reduced the 
"sagittal depth" of the curved mirror substrates and al- 
lowed us to maintain a gap of ^ 100 fj,m around the 
edge of the cavity. Machining of the mirror substrates 
was performed after they had been superpolished and 
coated, but this process did not seem to degrade the 



mirror reflectivities significantly. The nominal combined 
transmission and loss per mirror, before machining, was 
1.5 X 10~^. Direct measurements of the cavity finesse 
yield T ~ 217, 000 at an optical wavelength of 852.36 
nm, consistent with a combined mirror transmission and 
loss T* ~ 1.45 X 10^^. This value of T* is inferred from 
the measured value of I and the measured cavity HWHM 
k/2t: ~ 3.21 MHz. 

The cavity used throughout the work described in this 
paper suffers from a rather pronounced birefringence, 
which for TEMoo modes near 852 nm induced a splitting 
of 8±2 MHz between linearly-polarized eigenmodes. It is 
not entirely clear whether this birefrigence is a result of 
the substrate machining for diameter-reduction, a prop- 
erty of the coatings, or something associated with the 
mirror-mounting procedure used for this particular cav- 
ity. We note that recent efforts by other members of our 
group pl| ] have produced a cavity of approximately the 
same finesse with greatly-reduced birefringence (by a fac- 
tor ^ 25), using mirrors from a different coating run and 
with great care taken to minimize cavity misalignments 
and stress on the mirror substrates. 

Our two mirror substrates are mounted in vee-grooves 
atop independent aluminum blocks, with a piezoelectric 
actuator between the blocks for active servo-control of 
the mirror separation (cavity length). The cavity mount 
sits on a stack of alternating OFHC copper blocks and 
viton o-rings for passive vibration isolation, all within an 
ion-pumped vacuum chamber whose background pres- 
sure was typically ~ 10~^ torr (inferred from the ion 
pump current). In order to bring the MOT as close as 
possible to the central axis of the physics cavity, we had 
to use a rather open (and non-magnetic) mount design, 
leading to some compromises in the way of mechanical 
stability. Sitting on the vibration-isolation stack and un- 
der vacuum, we found that the native noise spectrum of 
the cavity length extended out to about 4 kHz, with one 
prominent resonance at 50 Hz (which we attribute to a 
transmission resonance of the isolation stack). Due to a 
set of PZT-actuator resonances above 10 kHz, we have 
ultimately been limited to a unity-gain bandwidth ~ 1 
kHz for the cavity stabilization servo. 

B. Laser and cavity locking schemes 

Although the principle aim for this experiment is to 
stabilize the cavity length at some precise offset A/27r ^ 
- 100 MHz from the Cs D2 resonance at 852.359 nm, 
the strong atom-cavity coupling places severe restrictions 
on the optical power that can be used for the purpose of 
generating an error signal. On resonance, the saturation 
intracavity photon number for our cavity is as small as 
mo = 75^/2^0 ~ 0.1, which sets a fiducial cavity through- 
put of 47rKbmo ~ 1 pW. With such low optical power it 
would be extremely difficult to obtain a high-quality error 
signal for locking the physics cavity. Other experiments 
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performed in our group have circumvented this problem 
by using a chopped locking scheme, in which a strong 
"lock beam" alternates with a weak "probe beam" at 50% 
duty cycle and ~ 1 - 4 kHz frequency [||,|o|j3|]. How- 
ever, such a strategy inherently limits the servo unity- 
gain bandwidths to ~ 100 Hz at best, and would not be 
appropriate for future experiments with atoms trapped 
inside the cavity for long periods of time. 

In this work we have developed an alternative locking 
scheme for the physics cavity, which employs an auxil- 
liary diode laser at 836 nm to monitor the cavity length 
on a different longitudinal mode than that which couples 
strongly to intracavity Cs atoms. We use a commercial 
grating-stabilized diode produced by New Focus (Santa 
Clara, California). With a 16 nm detuning, we can send 
~ 40 nW through the cavity and incur an AC Stark shift 
of only ~ 60 kHz for the atomic resonance at i^a- Using 
an EG&G avalanche-photodiode/transimpedance ampli- 
fier module (model C30998) for AC detection of transmit- 
ted 836 nm light, we obtain an FM error signal (modula- 
tion frequency 3.8 MHz) with signal-to-noise ratio (SNR) 
~ 20 at 30 kHz bandwidth. 

The 836 nm diode laser is stabilized to an auxilliary 
"transfer cavity," which consists of a pair of 25 cm radius- 
of-curvature mirrors at ~ 16 cm separation. One of the 
mirrors is mounted on a piezoelectric actuator to allow 
cancellation of DC drift and low-frequency noise. The 
transfer cavity has a linewidth ~ 100 kHz at both 836 nm 
and 852 nm, with an overall mode spacing ^ 300 MHz. 
The transfer cavity is also used for pre-stabilization of 
the Ti:Sapphire laser, and some of the TiiSapphire light 
is used in a Cs modulation-transfer spectrometer ||39| ] 
to provide an absolute reference for the transfer cavity 
length. From run to run we used one or two acoustoop- 
tic modulators to offset the Ti:Sapphire/transfer-cavity 
lock point by 4-140, +87, or +43 MHz relative to the Cs 
6Si/2(F = 4) -> 6P3/2(F' = 5) transition. 

By comparison of the TiiSapphire and diode laser error 
signals in their respective locks to the transfer cavity, we 
infer that the relative rms jitter between them is < 10 
kHz. For both laser locks we use the Pound-Drever-Hall 
technique of detecting an FM signal in reflection from 
the tranfser cavity. The stability of the transfer cavity 
resonances with respect to atomic Cs lines in a vapor 
cell was such that we did not see any relative jitter be- 
yond the measurement noise in our modulation-transfer 
spectrometer (SNR 50 : 1 in 30 kHz bandwidth). 

The TiiSapphire stabilization employs two feedback 
loops, one with ^ 10 kHz bandwidth to the tweeter inside 
the ring laser and another with ^ 100 kHz bandwidth 
to a VCO-driven, double-passed acoustooptic modulator 
(AOM) just outside the laser cavity |^ , ^ . The mean 
frequency of the AOM is 76 MHz, and the error signal 
going to the VCO has a lower cutoff of ~ 10 kHz to pre- 
vent DC drifts. We note that use of the AOM is crucial 
for achieving high stability of the TiiSapphire frequency. 
The diode laser servo utilizes both feedback to the grat- 
ing PZT and direct modulation of the injection current. 



achieving an overall unity-gain bandwidth 1 MHz. 

Having locked both the diode laser and TiiSapphire 
to modes of the transfer cavity, which itself is locked to 
Cs, we use a travelling-wave electrooptic modulator to 
generate an rf sideband of the diode laser at /o ~ 200 — 
500 MHz. Either the upper or lower sideband is used to 
derive an FM error signal for locking the physics cavity 
by dithering /o at 3.8 MHz, thus allowing us to achieve 
arbitrary placement of the physics cavity mode near 852 
nm via the tunability of /q. 

C. Evaluation of the physics-cavity servo 

Our basic requirement for the quality of the physics- 
cavity servo was that relative jitter of the cavity res- 
onance and the probe laser frequency should not con- 
tribute a significant amount of noise in the heterodyne 
photocurrent. Hence the relevant comparison to make is 
between the noise in both quadratures of a demodulated 
beatnote and the photocurrent fluctuations produced by 
the local oscillator alone. At 200 kHz bandwidth and 
with a probe beam strength such that m= |(a)p = 1.5 
in the empty cavity, the standard deviations of the phase 
and amplitude quadratures of the transmitted probe 
beam were measured to be 1.01 and 1.39 (respectively) 
relative to those of the local oscillator alone. Note that we 
have estimated the standard deviation of the quadrature- 
amplitude signals produced by our optical local oscillator 
to be only a factor of 1.05 above the theoretical shot-noise 
limit (see below). We therefore believe that our overall 
excess noise factor (3 ~ 1.5. 

By taking some simultaneous recordings of the hetero- 
dyne photocurrent and the physics-cavity error signal, we 
were able to verify directly that the atomic transits do 
not affect the physics cavity servo. 

D. Cesium MOT 

To provide a source of cold Cs atoms, we used a stan- 
dard magnetooptic trap loaded directly from a thermal 
beam BS. Our choice of thermal-beam loading, as op- 
posed to loading from a background vapor, was driven by 
an attempt to prevent accidental coating of the physics 
cavity mirrors with Cs. In more than two years of ser- 
vice, we did not detect any significant (> 5%) change in 
the cavity finesse. Pre-cooling of the Cs beam was not 
necessary for this experiment, as we required only a very 
low rate of delivering single cold atoms into the cavity 
mode volume. Our MOT employs a six-beam configura- 
tion, and we orient the anti-Helmholtz coils for the trap 
so that their symmetry axis is parallel to that of the op- 
tical cavity. This leads to a MOT laser beam geometry 
with one beam axis running parallel to and just above 
the cavity, plus two beam axes in the plane of the mirror 
surfaces (Figure ||). 
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The light for the MOT was provided by a pair of 
grating-stabilized diode lasers (SDL 5421-G2), one tuned 
to the Cs 6Si/2(F = 4) ^6P3/2(F' = 5) cycling 
transition for trapping and the other to 6Si/2(^ — 
3) — S'6P3/2(^"' = 4) for repumping. Each trapping beam 
had ~ 1 cm diameter and anywhere from 40 /iW to 4 
mW of optical power, depending on how many atoms we 
were trying to send into the physics cavity. We typically 
used a Cs reservoir temperature of 60-80 C for the ther- 
mal beam, which effused through a 200 n pinhole and 
travelled an overall distance of ~ 60 cm to the trapping 
region (with a cold mechanical collimator in the way to 
reduce loading of the ion pump). With an anti-Hclmholtz 
field of around 25 G/cm, we could load up to ^ 2 x 10^ 
atoms into a millimeter-sized cloud, whose mean temper- 
ature we estimate to be ~ 100 /iK based on fluorescent 
imaging of free expansion. This temperature estimate is 
also supported by the spread in arrival times of individual 
atoms falling into the cavity. 

When running the experiment we would load the MOT 
for about 0.5 s, then drop it by quickly turning off 
the trapping beams with an AOM (using an rf switch 
with 45 dB attenuation). After the trapping beams 
were thus extinguished, we would ramp down the anti- 
Helmholtz field according to an RC-filtered step with ~ 4 
ms time constant. The repumping beam was left on all 
the time, so that falling atoms would be shelved in the 
F = A ground hyperfine level before entering the cavity. 
No specific preparation was performed with respect to 
the atomic Zeeman states. Dropping ~ 2 x 10^ atoms we 
would generally see 30-50 atoms falling through the cen- 
tral part of the cavity mode volume, so for single-atom 
transit data we had to reduce the trapping beam power 
by a factor of 100 to reach 0-2 atoms per drop. The over- 
all repitition rate for the trap-drop cycle was typically 
0.6 Hz. 



E. Probe generation and photodetection 

We used a balanced-heterodyne setup in order to 
achieve high-eflSciency, zero-background photodetection 
of ~ 1 pW levels of 852 nm light transmitted through the 
physics cavity. The frequency difference between cavity 
probe light and the optical local oscillator for heterodyne 
detection was between 40-190 MHz, depending upon our 
choice for the atom-probe detuning. The probe light was 
generated from the Ti:Sapphire output by cascading a 
-1-200 MHz AOM and a tunable travelling- wave electroop- 
tic modulator, which was driven between -245 MHz and 
-440 MHz to produce the desired atom-probe detuning. 
This indirect method was required to prevent contamina- 
tion of the heterodyne photocurrent by electronic noise 
at the heterodyne frequency. 

Light leaving the physics cavity first hit a color- 
separation mirror which reflected > 99% of the 852 nm 
light but transmited ~ 30% of the 836 nm light, allowing 



us to recover an error signal for locking the physics cavity 
(see above) without compromising the overall detection 
efficiency for the probe field. Residual 836 nm light going 
to the heterodyne setup amounted to only ~ 30 nW and 
had negligible effect on the photocurrent of interest. 

The local oscillator (LO) for the optical heterodyne 
was spatially cleaned by a J-" ~ 1000 Fabry-Perot cavity 
(linewidth ~ 1 MHz) , which also served to strip off spec- 
tral noise at 76 MHz associated with the AOM servo for 
stabilization of the Ti:Sapphire frequency. The cleaning 
cavity was locked using the Pound-Drever-Hall method 
[0 with FM sidebands at 24 MHz, which likewise had 
to be kept weak in order not to saturate the AC gain of 
the heterodyne photodetectors. We used a total of ^ 2 
mW in the LO, which generated a shot-noise level ^ 5 
dB above the electronic noise of the photodetectors in 
the frequency range of interest. 

The difference photocurrent from the balanced hetero- 
dyne detectors was amplified up to -50 dBm or higher, 
then divided by a 0° rf splitter. An independent signal 
generator was used to produce an rf local oscillator at 
the heterodyne frequency, and it was halved using a 90° 
rf splitter. The two identical copies of the photocurrent 
were mixed with the in-phase and quadrature copies of 
the rf LO to produce an orthogonal pair of quadrature 
amplitude (QA) signals at baseband. The QA signals 
were further amplified, and passed through 300 kHz ana- 
log filters with a roll-off of 12 dB/octave. We used a 
12-bit ADC to sample both QA's simultaneously at a 
rate of 10 MHz per channel, which is sufficiently high to 
avoid signal aliasing completely. Following each drop of 
the trap, we continuously recorded both QA's for a data 
acquisition window of 50 ms and streamed the data to a 
hard drive for offline processing following the experimen- 
tal run. 

Ideally, we would like the data acquisition procedure 
just described to yield directly the amplitude and phase 
quadrature amplitudes of light transmitted through the 
cavity. If we write the transmitted optical fleld as 
^(i)e^"^p*, where A{t) is a slowly- varying complex am- 
plitude, the amplitude Xa and phase Xp quadrature- 
amplitudes are defined hy A= Xa + ixp. With respect to 
the quantum-mechanical theory of the Master Equation 
(|l|) , A (X {a). Note that we define A to have zero phase 
when the cavity is empty, so that Xp should have zero 
mean when there are no intracavity atoms. 

Given the way that we generate the probe beam, how- 
ever, we have no way of generating a phase- locked rf local 
oscillator to recover Xa and Xp directly. The phase of the 
heterodyne photocurrent differs from the phase that the 
light has just after it leaves the cavity because of fluctu- 
ations in the relative optical path length travelled by the 
signal beam and optical local oscilator in reaching the 
photodetectors. So the two signals produced by mixing 
the photocurrent with the shifted and unshifted copies 
of our rf local oscillator correspond to an orthognal, but 
rotated pair of quadrature amplitudes xi, X2' 
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xi \ ^ f cos(j> -sin(l)\ I' Xa\ 
\X2 J ysmcf) cos 4> J \Xp J ' ^ ' 

Luckily, the characteristic timescales for fluctuations in 
the phase (f> are quite long (~ 10— 100 ms, corresponding 
to acoustic disturbances) compared to the 250 fis dura- 
tion of an individual atom-transit signal. In processing 
the recorded data to produce the plots discussed below, 
we have therefore used an "adaptive" definition of the 
amplitude and phase quadrature amplitudes. Within a 
window of 2 ms preceding the signal of interest, we es- 
timate the instantaneous value of (j) by determining the 
rotation of xi,X2 that produces one quadrature Xp with 
zero mean and one quadrature Xa with positive mean. 
Then Xp is operationally defined to be the phase quadra- 
ture photocurrent, and Xa is the amplitude quadrature 
photocurrent. 

Figure shows an example of a 15 ms segment of our 
quadrature amplitude data. Note that some excess low- 
frequency noise can still be seen in the phase quadrature, 
which constrains the lower end of our measurement band- 
width to ~ 1 kHz. Six prominent atom-transit signals, 
characterized by a sharp drop in the amplitude quadra- 
ture and a simultaneous increase in the phase quadrature, 
can be seen between t — 0.009 sec and t = 0.013 sec. In 
our subsequent discussions of the data, we shall focus on 
individual signal "events" of this type. 

III. DETERMINATION OF INTRACAVITY 
PHOTON NUMBER 

In order to make quantitative comparisons between 
data and theory, we need to calibrate the strength of the 
driving field used in each experimental run. As we will 
ultimately choose to ignore DC optical phase offsets, the 
relevant quantity for us will be the number of photons 
m that builds up when the cavity is empty. This corre- 
sponds to 2Ka\£\'^ / (na + Kft + Kcj in the Master Equa- 
tion case and \y\^ /mo for the OBSE. 

Our strategy for determining the intracavity photon 
number during experimental runs has been to work back- 
wards from the heterodyne photocurrent observed with 
no atoms in the cavity, using the expression ||3^ 

(^)^_^_=4T,..™. (10, 

where S is the mean value of the demodulated amplitude 
quadrature signal (i.e. the heterodyne photocurrent is of 
the form S cos{ujt + (j))), N is the mean-squared power of 
fluctuations in the same signal due to (optical) shot noise, 
T is the measurement interval, and m = | (a) p. Note that 
the above expression is valid for a coherent signal beam, 
which we assume to be the state of the light transmitted 
through the cavity when no atoms are present. For an 
accurate calibration, we thus need to know the output 
mirror transmission Kf,, the total cavity loss (which can be 



determined from measurements of the cavity finesse), and 
the overall photodetection efficiency 77. In determining m 
from the data, we typically chose T on the order of 1 ms. 
The following subsections provide further detail for each 
critical aspect of the calibration. 

A. Evaluation of local oscillator noise 

Ideally, the quantity N appearing in equation ( |l0| ) 
should coincide with the photocurrent noise power n ob- 
served when the signal beam is blocked. This would allow 
us to determine the intracavity photon number without 
having to calibrate the exact gains of the photocurrent 
amplifiers, etc. But the procedure is invalid if the opti- 
cal local oscillator carries excessive technical noise, so we 
have empirically checked the scaling of our nominal N 
with the DC optical power of the LO. A fit of the data 
to the functional form n — aP + bP^ yields b/a ~ 0.11 
mW^^. Given our typical operating LO power of 1 mW 
per detector, and considering the relative magnitudes of 
other uncertainties, we approximate TV ~ n. Note that 
the LO does pass through an intensity stabilizer with 
400 Hz bandwidth just before entering the cleaning 
cavity. 

B. Measurement of cavity decay rates 

We directly determined the total cavity field decay 
rate Ka + Kb + Kc by making a calibrated measurement 
of the HWHM of a cavity TEMqo resonance. In order 
to minimize systematic errors, we did this by using the 
heterodyne detectors to monitor the transmitted optical 
power of a fixed-frequency Ti:Sapph probe beam while 
scanning the cavity length. The cavity length was al- 
ways under servocontrol during the measurement, as we 
generated the scan by stepping the frequency of the rf 
going to the travelling wave modulator for the locking 
diode laser. The total cavity field decay rate is then 
given by the resonance HWHM measured in terms of the 
modulator rf scan, times a correction factor of the ratio 
of the diode laser and Ti:Sapph wavelengths. We find 

Ka ~t- Kb -\- Kq — 3.2 MHz. 

As the mirrors used to construct the cavity should be 
identical, we assign half the total losses to each mir- 
ror. Unfortunately, we did not manage to characterize 
the ratio of intracavity losses Kc to transmission losses 
Ka + Kb before we accidentaly damaged the cavity. Given 
that the mirror coatings we have previously received from 
the same manufacturer have displayed very low scatter- 
ing/absorption loss (~ 1.1 X 10~^) in the wavelength 
range of interest |Q, we have assumed Kb — 1.6 MHz. 



6 



C. Measurement of rj 

Three principal factors determine the value oirj = Vet: 
the spatial overlap -s/V between the signal beam and the 
optical local oscillator, the photodetector quantum ef- 
ficiency e, and the (power) efficiency r with which we 
transfer light from the cavity ouput to the photodetec- 
tors. Using an optical powermeter we measured r — 0.9, 
as well as £ ~ (0.68, 0.71) for our two photodetectors. 
We measured V ~ 0.65 by producing a DC fringe be- 
tween the signal beam and the optical local oscillator, 
with both beams adjusted to a power level of 1 /iW. The 
power adjustments were made upstream of the physics 
cavity and LO cleaning cavity, ensuring that no optical 
misalignments were introduced in the process. We ob- 
tained an independent measurement of Ve ~ 0.4 ± 0.05 
from the heterodyne signal-to-noise ratio obtained with 
a measured signal power of 10 nW. Note that the reason- 
able agreement between this latter value and the directly- 
measured Ve provides further indication that our optical 
local oscillator bears minimal excess noise, and that our 
assumptions Kb = Kc ~ are valid. 

Although it could have been avoided, we did suffer one 
additional loss of detection efficiency due to gain and 
phase imbalance in the subtraction of heterodyne pho- 
tocurrents. Several weeks after we took the data sets 
shown below, we realized that one photodetector output 
had an rf signal 4 dB higher than the other (this figure 
includes both the deviation of our heterodyne beamsplit- 
ter from being 50/50, the difference in photodiode quan- 
tum efficiencies, and the difference in transimpedance 
gains), and a phase offset of 1 radian. If we write 
the two photocurrents (for a coherent signal beam) as 
ii = Tje-^"^* + a and ^2 - ^6"*'^*+*'^ + (where Ci 
and a are uncorrelated gaussian white noises with zero 
mean and unit variance), the difference photocurrent is 
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where ^ is again a gaussian white noise process with zero 
mean and unit variance. The complex imbalance ^e*"^ 
should thus reduce the effective photodetection efficiency 
(for the purpose of evaluating the heterodyne signal-to- 
noise ratio (p^) by 
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Using the measured g ~ 0.63 and (j> 
and we quote an overall value of 77 : 



(12) 

: 0.85 rad, 77 0.8??, 
0.32. 



atoms fall through the cavity.]^ Our code simulates three- 
dimensional, classical center-of-mass motion of individual 
atoms under the influence of fluctuating forces due to 
their strong interactions with the cavity field. Gravity is 
also included in the kinematics, but plays only a minor 
role over the ~ 1 ms duration of the simulations. 



A. Overall scheme for the simulations 

The basic scheme of the simulations is to precompute 
the values of the mean cavity-induced force, the cavity- 
field-induced diffusion coefficient, and the steady-state 
(complex) amplitude for the intracavity field as a func- 
tion of atomic position. Hence, every other degree of 
freedom in the simulation is adiabatically eliminated and 
slaved to the atomic motion, which is assumed to be the 
slowest and the "stiffest" process in the dynamics. The 
mean force and intracavity field are derived from steady- 
state solutions of the Master Equation (^, and the cav- 
ity diffusion is computed using the Quantum Regression 
Theorem |Q. In each timestep, the code first performs 
an interpolation on the precomputed tables of values to 
determine the appropriate change in atomic momentum 
(which includes a stochastic increment consistent with 
the local value of the diffusion constant), records an ap- 
propriate value for the cavity output field, and then up- 
dates the atomic position and momentum. The simplest 
possible integration scheme is used (explicit Ito-Euler), 
yielding order 0.5 convergence in the timestep. We chose 
a timestep of 7.5 ns in order to keep the run-times for the 
simulations reasonable under Matlab on our Pentium II 
workstations, and this should have been sufficiently small 
to keep the integration error below other sources of inac- 
curacy. 

The principle shortcoming of this scheme is that it as- 
sumes the atomic velocity will remain small enough that 
variations in the coupling strength g will be negligible 
over timescales on the order of and 7"^ . If this condi- 
tion is violated, then the steady-state values of the mean 
cavity force and cavity output are no longer appropri- 
ate. Another major approximation was made in treating 
only a two- level atom — optical pumping and/or the asso- 
ciated opto-mechanical effects {e.g., spontaneous forces) 
could certainly play an important role in determining the 
shapes of our observed transit signals. A third, less severe 
approximation is that we have treated the stochastic pro- 
cess associated with recoil from spontaneous-emissions as 
being statistically independent from the one associated 
with dipole-force fluctuations. The possibility that in- 
teresting effects could arise due to quantization of the 
atomic motion is of course ignored as well. 



IV. NUMERICAL SIMULATIONS 



In order to facilitate the interpretation of our data, 
we generated some rudimentary Monte Carlo simula- 
tions of the heterodyne signals that we should see as 



^The computational work described in this section was con- 
ducted by H.M.. 
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Given the relative simplicity of our simulation scheme, 
it would be inappropriate to draw any strong conclu- 
sions about our experimental data soley on the basis of 
apparent similarities with the numerical results. How- 
ever, the simulations do provide an appealing model for 
the atomic center-of-mass dynamics and associated het- 
erodyne signals, with predictions that seem to be fully 
consistent with what we see in the actual experiment. 
This model suggests some interesting interpretations for 
qualitative features of the data. For example, the marked 
asymmetry (in time) of most of our experimental transit 
signals may be associated with an initial phase of atomic 
"channeling" in the dipole potential provided by the in- 
tracavity optical field, followed by a sudden escape from 
local confinement due to cavity-field-induced momentum 
diffusion. Such a dynamical process, if it could be con- 
firmed in some way, would bear a strongresemblence to 
the Levy Walk behavior predicted in and observed 

in 1^^. Our simulations also suggest that isolated dips 
and steps in the observed transit signals may be associ- 
ated with sudden changes in atomic localization relative 
to the cavity standing wave. 

Before presenting some results of the atomic-motion 
simulations, let us first describe the numerical method 
used to compute the cavity diffusion coefficients. 

B. Computation of diffusion coefficients 

Following Doherty et al we computed the diffu- 
sion coefficient associated with dipolc-force fluctuations 
according to 

D= lira Re {{{t),i{t - T))dT, (13) 

where f is the force operator 

f = -ifiWg (r) (at6-_ - da+) (14) 

and we define 

g (r) = 50 cos{kLx) exp [- {y^ + z^) /w'^] . (15) 

Note that the vector nature of the force operator comes 
only from the gradient (r). Unlike the computation of 
mean forces and the expected intracavity field amplitude, 
the evaluation of (|l3|) requires an actual time-integration 
of the cavity-QED Master Equation (|l|). Note that 

(f(t),f(t-r)) ^Tr[fe^-fp,,] -(f)2, (16) 
where C is the Liouville superoperator defined by 

j^P = Cp, (17) 

with reference to the cavity-QED Master Equation (|^). 
Knowing the steady-state atom-cavity density matrix 
Pss, one can evaluate the correlation function in (|l^) by 



integrating the Master Equation for a time r with fpss 
as an initial condition, multiplying the result by pss, and 
finally taking the trace. 

The numerical integration was performed by a 
C-H-h/MPI code on an SGI/Cray Origin-2000 cluster. We 
used a truncated basis of 25 Fock states for the cavity 
mode, using an explicit Euler integration with a 1 ps 
timestep and 5 ps total integration time. 

Figure ^ shows the results of such a calculation for the 
diffusion coefficient as a function of atomic position along 
the standing wave, for an atom located on the cavity axis. 
Recall that the cavity mode function varies as cos{kLx) 
along the cavity axis, and as a gaussian in the transverse 
dimensions. Starting from the top, the four curves repre- 
sent A = 0, 10, 30, and 50 MHz, all with a probe strength 
such that an average of 2 photons would build up in the 
cavity if it were empty and with = 0. Figure |5| shows 
similar results, all for A = 50 MHz, but with an aver- 
age intracavity photon number (for the empty cavity) of 
m = 8 for the top curve, 4 for the middle, and 2 for the 
lowest. As the variation of g along the standing wave is 
much more rapid than in the transverse (Gaussian) di- 
rections, we approximated V.g (r) ~ dg/dx in expression 
(p^). The computation of f for the mean force field was 
of course three-dimensional. 

Note that expression (^3|) gives the diffusion coefficient 
associated with dipole force fluctuations only — it takes no 
account of recoils from spontaneous emission. We there- 
fore incorporate a second diffusion process in the simula- 
tions, whose coefficient is computed from the expectation 
value of the atomic excitation at each point in space: 

Drec=^^^r{a+a^). (18) 

Here T is the atomic spontaneous emission rate (Einstein 
A coefficient), and the number 1/25 comes from averag- 
ing over the angular distribution pattern for dipole radi- 
ation ]50| ]. 

C. Discussion of results 

Figure || shows an example of the atomic position 
along the standing-wave x and the atomic velocity along 
the standing-wave from one typical simulation. Also 
shown are the mean value of the intracavity phase 
quadrature amplitude (qp) = i{a — a^)/2, and a sim- 
ulated heterodyne signal with Gaussian noise added at 
the level appropriate to our experimental parameters (see 
equation (|lO|)). Note that the quantity {{qp) plus noise) 
is proportional to the photocurrent Xp that would be 
recorded in our experiment. In the bottom row of the 
figure, {qp) and {qp) plus shot-noise are shown at an ana- 
log bandwidth 300 kHz in order to illustrate the effects 
of finite detection bandwidth on the qualitative features 
of the data [Q. The results clearly suggests that the 
gross features of the simulated signals are most strongly 
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affected by atomic motion along the standing wave, in 
that the overall gaussian profile associated with motion 
in the z-direction becomes distorted by the "envelope" 
of the oscillatory variations due to motion along x (recall 
g (r) = cos{kLx) exp [— (y^ + z^) /w^] ). 

Looking at the results for x{t) and Vx{t), we see that 
the atom in this particular simulation was mechanically 
confined within one well of the cavity standing wave until 
t ^ 810 /xs. That is, x did not vary by more than A/4 
and Vx displays the oscillations that one would expect to 
see for an atom trapped within a potential well. At time 
t ~ 810 ^s, however, we see that momentum diffusion fi- 
nally pushes Vx above some threshold value such that the 
atom is able to escape from local confinement and "fly 
above" the periodic dipole potential associated with the 
cavity standing wave. Generally speaking, one expects 
this type of escape to occur as the atom passes through 
the cavity axis in the — z-direction and the standing- wave 
potential wells begin to decrease in depth. The expec- 
tation value {qp) clearly reflects the qualitative features 
of the atomic trajectory x{t), although we also see that 
the details are lost in the measured signal due to het- 
erodyne shot-noise. What survives in the simulated het- 
erodyne signal is an overall asymmetry in the envelope 
of the transit-signal, with a sudden "step" downwards at 
t ^ 800 /xs (as suggested by the dotted vertical line). 

Figure ^ shows four simulated signals with A = 10 
MHz and an average intracavity photon number of m = 
2, and two simulated signals with A = 50 MHz and 
m = 4 (6 = in all cases). In each subplot, the up- 
per trace shows the intracavity amplitude quadrature 
{qa) = {a-\- a)) /2 versus time, and the lower trace (with 
zero mean value) shows the intracavity phase quadrature 
{qa) versus time. Again, Gaussian noise has been added 
to the simulated signals in order to reproduce the overall 
signal-to-noise ratio predicted by equation (p^. In all 
simulations, the atom was started at a position 7 gaus- 
sian waists above the cavity with a vertical velocity of 
-47 cm/s. The initial transverse position, transverse ve- 
locity, position along the cavity standing-wave axis, and 
velocity along the standing- wave were randomly assigned 
within parameter ranges that were narrow enough to en- 
sure that most simulated atom-drops produced a sizeable 
signal. 

The simulated signals in Figure |^ may be compared 
directly to the real data shown in Figure ^. Like the sim- 
ulation shown in Figure the signals in Figure ^ display 
a generic asymmetry and often contain an isolated dip 
or step. The features of this type that are marked by an 
arrow near the horizontal axis are again associated with 
sudden changes in the atomic localization relative to the 
cavity standing wave (as determined by examining x{t) 
from the simulations). 

Although we have the beneflt of knowing both the "ac- 
tual" trajectory of an atom and the corresponding het- 
erodyne signal in our simulations, a priori knowledge of 
the atomic trajectory is of course unavailable in our ex- 
periment. Hence, the results of our numerical simulations 



can only provide general guidelines for how we might try 
to infer something about individual atomic trajectories 
from our heterodyne transit signals. First of all, it ap- 
pears that signflcant asymmetries in the observed signals 
could provide evidence that we do in fact observe me- 
chanical aspects of the atom-cavity coupling within in- 
dividual atomic transits. Relative to the approximations 
involved in our simulations however, one should bear in 
mind that such asymmetries could also arise from optical 
pumping processes that may occur during at atomic tran- 
sit through the cavity. Second, the simulations suggest 
that the transit signal envelopes are most strongly influ- 
enced by atomic motion along the cavity standing wave, 
as opposed to details of the motion along y oi z. Third, 
isolated dips or steps in the signal could be indicative of 
sudden changes in the atomic localization with respect 
to the standing wave. Again however, we should be care- 
ful to note that such features might also be caused by 
intracavity optical pumping. 

V. DATA 

Moving on to the experimental data, let us first dis- 
cuss some transit signals displayed in the same fashion 
as the simulations of Figure 0. The quality of this data 
illustrates our experimental ability to perform continu- 
ous, nearly quantum-noise limited measurement of the 
atom-cavity interaction energy during individual scatter- 
ing events. In the second subsection we shall display tran- 
sit data on the complex plane, and present a quantitative 
comparison to theoretical predictions of the atom-cavity 
system's complex optical susceptibility. We shall see that 
in optimal cases, sufficient information may be obtained 
from an individual atom-transit signal to distinguish un- 
ambiguously between quantum and semi-classical models 
of cavity QED. 

Note that all of the data shown here were taken with 
cavity-probe detuning Q = Vc — Vp ~ 0, while the atom- 
probe detuning = Va — and probe power m are 
varied. 



A. High bandwidth single-atom transits 

The basic unit of our experimental data is a continu- 
ous stream of quadrature-amplitude values versus time, 
recorded for ^ 50 ms windows following each dropping of 
the MOT. One ~ 15 ms segment of such data is shown in 
Figure ^. This particular data segment was taken with 
detunings (A = 10, O = 0) MHz and with a probe power 
such that an average of 1.5 photons would build up in 
the cavity if it were empty (m — 1.5). 

Figure || shows six of the largest single-atom transit 
signals from our entire data set, which covers detunings 
A from —10 to -1-100 MHz and probe powers m from 1.5 
to 11 photons. The particular values of A and m for each 
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signal in Figure || are displayed by the y-axis. The data 
are shown at our full analog bandwidth of 300 kHz and 
sampling rate of 10 MHz (12 bit resolution). Note that 
we have displaced the amplitude quadrature signals Xa 
by +400 in order to prevent them from overlapping with 
the phase quadrature signals Xp. 

We have found that transit signals of maximal contrast 
tend not to have very much internal structure, although 
there is a fairly prominant dip at the point indicated 
by an arrow in subplot "dat4," a set of three dips in- 
dicated by the arrow in subplot "dat6," and oscillatory 
structure in the signal of subplot "dat5." The signal 
of subplot "dat4" is really not so different in its overall 
structure than the simulations shown in Figure ^ and in 
the "trajlO" subplot of Figure |. One can see that the 
shapes of the overall signal envelopes do vary substan- 
tially and are generally asymmetric in time. The distri- 
bution of the atom-transit signal between amplitude and 
phase quadratures clearly depends upon the probe de- 
tuning, as will be discussed in greater detail below. For 
the data event recorded with A = 0, one sees that there 
is only a reduction in the power of the transmitted probe 
without any shift in its phase. At all other detunings, 
the atom-induced reduction of amplitude quadrature is 
partly offset by an increase in the phase quadrature, indi- 
cating a significant phase shift of the transmitted probe 
beam. In subplot "dat6" taken with A = 50 MHz, the 
signal is seen to primarly reside in the phase quadrature. 

For signals such as those shown in Figure ^ we may 
estimate a full-signal to rms-noise ratio of approximately 
2.5 for the phase quadrature and 4 for the amplitude 
quadrature. Combining the two signals, we have an over- 
all ratio of 4.5. The signal bandwidth is 300 kHz, imply- 
ing a relative sensitivity of 0.22/V3 x 10^ ~ 4 x lO"'' 
Hz~^/^. Assuming that the largest signals correspond to 
atoms reaching the maximal coupling strength of 50 = H 
MHz, this sets our broadband sensitivity Sg to time- 
varations in g = Ei^t/h to be Sg ~ 4.5 kHz/\/Hz. Given 
our overall photodetection efficiency 77 ~ 0.32 and our 
independent assesment of the excess noise factor /3 ~ 1.5 
(corresponding to technical noise on the optical local os- 
cillator power and on the physics cavity lock), we find 
that this value lies only a factor of P/^/rj ~ 2.7 above the 
fundamental quantum noise floor. We maintain this high 
sensitivity over 2.5 decades of signal frequency, from our 
full photodetection bandwidth of 300 kHz down to 1 kHz 
(at which point we are limited by residual technical noise 
in our measurements of the phase quadrature). 

Although a complete theoretical description of our 
experiment at the level of a Stochastic Master Equa- 
tion [ll4| could in principle be formulated to include 
three-dimensional quantized atomic center-of-mass mo- 
tion, cavity birefringence, and the full manifold of in- 
ternal atomic Zeeman states, the resulting model would 
almost certainly be too complex to be useful for quan- 
titative interpretations of our existing data. Neverthe- 
less, we wish here to stress that only a theory that ex- 



plicitly treats the quantum conditioning [|l^-|T^ of the 
atom-cavity evolution on the stochastic component of the 
heterodyne photocurrent could provide a full account of 
the statistics and autocorrelation properties of our atom- 
transit data. In this crucial sense, we would claim that 
our recent work may be distinguished from all other ef- 
forts in contemporary experimental physics. Of course, 
much remains to be done before we can offer definitive 
proof, so our continuing efforts will largely be directed to- 
wards the unambiguous demonstration of quantum con- 
ditional dynamics in cavity QED. 



B. Transit phasors 

In addition to displaying our atom-transit signals in 
the format of photocurrent versus time, we can also con- 
struct parametric plots of Xa versus Xp to examine the 
correlation induced between these two quantities by the 
atom-cavity interaction. This is equivalent to compiling 
a histogram of the complex amplitude of the optical field 
transmitted through the cavity during an atomic transit, 
and hence to a continuous monitoring of the atom-cavity 
system's complex optical susceptibility. The underlying 
probability distributions for such histograms are deter- 
mined by the (evolving) Husimi Q-function for the intra- 
cavity field 

Many of our measurements have been conducted in a 
regime of strong cavity-driving fields (m > mogo/^^)i 
with the consequence that significant effects of satura- 
tion and optical nonlinearity can be seen in our exper- 
imental data. As has been extensively discussed in the 
cavity-QED literature, the nonlinear optical response of 
the atom-cavity system represents an important exper- 
imental signature that may be used to distinguish be- 
tween quantum and semiclassical theoretical models for 
the atom-cavity interaction. Prior to this work, four 
groups [pT|-p3 31 1 have reported nonlinear measurements 
in cavity QED. In each of these previous experiments, 
the information gained per atomic transit about the de- 
tails of the system's nonlinear response was far less than 
in the data we shall present here. Hence, we wish again 
to stress that the techniques developed in our most re- 
cent experimental work have opened new possibilities for 
exploring quantum dynamics at the level of single quan- 
tum realizations (trajectories), as opposed to the level of 
ensemble averages. 

In Figures ^ and |l^, we display single-transit data on 
the complex plane. The experimental data are displayed 
as gray dots, with each dot representing the values of Xa 
and Xp at some particular time during an atomic transit. 
In order to produce a set of dots from a transit signal such 
as those shown in Figure ^, the quadrature amplitude 
signals were first passed through an anti-aliasing filter 
and then resampled every 10 /iS to produce a discrete set 
of points {xl^jXp }. For every i, a gray dot is placed at 

radius \ {x\^Y + {xpY and polar angle [Xp / x'^^) . 
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These transit phasors represent a simple way of look- 
ing at the correlation between the ampHtude and phase 
of light transmitted through the cavity, with time re- 
moved from the picture. We beheve that this type of 
plot is the best format for comparisons to theory, be- 
cause atomic internal (optical pumping among Zeeman 
states) and external (motion through the cavity eigen- 
mode) dynamics should be factored out. In the simplest 
approximation, these two types of processes merely in- 
duce fluctuations in the atom-cavity coupling g, which 
should only move the locus of our gray dots in and out 
along a curve parametrized by g. The overall shape of 
this curve should be dictated by the interaction Hamilto- 
nian for the atom and cavity mode, and is therefore quite 
easy to compute. 

In Figure O, we show a set of transit phasors taken at 
various probe detunings. The data overlay theoretical 
curves predicted by the quantum and semiclassical theo- 
ries. The quantum-mechanical predictions are computed 
by finding the steady-state solution of the Master Equa- 
tion (P for values of g in the range [0, go] and appropriate 
values of A and m. The solid curves represent an interpo- 
lation through the discrete set of computed values for (qa) 
and {qp) as a function of g. The semiclassical prediction 
is computed in the analogous manner, using the Optical 
Bistability State Equation rather than the quantum 
Master Equation. Note that both curves must agree for 
g — 0, and this common point of origin is marked in each 
subplot by a triangle. The g = go endpoint of the quan- 
tum curve is marked by a circle, and the g = go endpoint 
of the semiclassical curve is marked by an x . In the sub- 
plots of Figure ^, the grey dots represent an overlay of 
data from two individual atom-transit signals. The two 
particular signals displayed in each subplot were selected 
on the basis of having maximal "contrast" among all the 
signals from a given data set, under the assumption that 
the atoms causing maximal signals should have come the 
closest to actually achieving g = go on their way through 
the cavity. Quantum theory is seen to predict the obser- 
vations quite well for all the detunings shown in Figure 

I- 

The quantum and semiclassical predictions shown in 
Figure || do not differ significantly except for the case of 
A = 0. Generally speaking, the two theories are known 
to agree in their predictions for weak-field response of 
the atom-cavity system but to differ in their predictions 
regarding saturation of the optical response. We here 
present experimental evidence from single atom transits 
for the quantum character of saturation in the atom- 
cavity system's response to near-resonant driving fields. 

In Figure ^ we show a set of transit phasors for fixed 
detuning A = 10 but variable probe strength m = 2.8, 
4.4, 7, and 11. A clear discrepancy between the quantum 
and semiclassical theory may be seen in this sequence 
of plots, with the data (gray dots) showing significant 
preference for the quantum-mechanical predictions. Each 
subplot in Figure |l^ represents an overlay of data from 
three individual atom-transit signals, in order to fill in 



the overall "shape" of the experimental transit phasors. 
It should be clear, however, that just a single data trace 
would suffice to determine that the experiment matches 
much more closely with the quantum theory than the 
semiclassical theory. We plan to elaborate this result in 
separate publication js^ . 

VI. CONCLUSIONS AND FUTURE GOALS 

In summary, we have described the details of our recent 
experimental work to perform continuous measurement 
of the interaction energy _Eint = Jig between one atom 
and an optical cavity during individual transit events. 
We displayed heterodyne transit signals in two comple- 
mentary formats, one of which highlights the large band- 
width and nearly quantum-noise limited signal-to-noise 
ratio achieved in tracking the time-evolution of the am- 
plitude and phase quadratures (ia,ip) of light transmit- 
ted through the cavity. The second data format (transit 
phasors) displays the correlation induced between Xa and 
Xp by the atom-cavity interaction, and our experimental 
results show that we are able to distinguish unambigu- 
ously between quantum and semiclassical models of cav- 
ity QED in only a few (or even just one) atomic transits. 
We have furthermore presented the results of rudimen- 
tary numerical simulations of atomic motion under the 
influence of mechanical forces and momentum diffusion 
associated with the strong atom-cavity coupling, and we 
examined the interpretations suggested by these simula- 
tions for certain qualitative features of our experimental 
data {i.e. steps and asymmetries). 

The primary conclusion we wish to draw in compar- 
ing the simulation results with our experimental data is 
that we have experimentally reached a regime of mea- 
surement sensitivity and bandwidth in which details of 
the atomic center-of-mass trajectories really should be 
visible. Even though we have no incontrovertible means 
of proving that steps and asymmetries seen in our exper- 
imental data should be associated with dynamical pro- 
cesses like channeling and diffusive escape, we are now 
motivated in our continuing work to develop some means 
of actively influencing the atomic motion while it is still 
inside the cavity. This would allow us to produce delib- 
erate displacements of an atom along the cavity stand- 
ing wave, and to examine the induced variations of the 
heterodyne signal in order to verify our inferred displace- 
ment sensitivity of 1.5 x 10~^° m/^/Hz. We anticipate 
that the combined abilities of monitoring and influencing 
atomic position relative to the cavity standing wave will 
enable the investigation of schemes for real-time feedback 
control of quantized atomic ccnter-of-mass motion. 

Finally, let us note that the standard quantum limit 
for overall observation time of the position of a single 
atom with sensitivity 1.5 x 10"^" m/\/Hz and 300 kHz 
bandwidth should be t, ~ 10 fis [|5|. This implies that 
quantitative experimental investigations of conditional 
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quantum dynamics (as described in pq]) should indeed 
become possible in our cavity QED system once we are 
able to reliably prepare individual atoms in well-defined 
initial states of motion. The ideal situation in this regard 
would be to trap and localize atoms within the cavity, re- 
leasing them at a node or antinode of the standing wave, 
on the cavity axis, and with an initial position uncer- 
tainty that is small compared to A/4. Current efforts in 
our group focus on trying to achieve this level of control 
via optical dipole- force traps and/or far-detuned optical 
lattices inside the cavity. 
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FIG. 1. Schematic overview of the apparatus. 

FIG. 2. Geometrical arrangement of the MOT beams rel- 
ative to the physics cavity (figure not to scale). The mirror 
substrates are each 3 mm in diameter and 4 mm long. The 
MOT forms at a height of ~ 7 mm above the cavity axis. 

FIG. 3. A 15 ms segment of the typical recorded data, 
showing the simultaneously-recorded photocurrents propor- 
tional to the amplitude quadrature Xa (upper trace) and phase 
quadrature Xp (lower trace) of light transmitted through the 
physics cavity (see text). The experimental parameters for 
this data were A = 10 MHz, 6 = 0, and m = 1.5. The 
photocurrents were digitized at a rate of 10 MHz, with 12-bit 
resolution, and the analog bandwidth of the anti-aliasing fil- 
ters was 300 kHz. Note the set of transient features clustered 
between 0.009 and 0.013 on the time axis, each of which was 
caused by the passage of an individual atom through the cav- 
ity. 

FIG. 4. Cavity diffusion coefficients versus atomic position, 
for (from the top curve going down) atom-probe detuning 
A = Ua — Vp = 0, 10, 30, and 50 MHz, with m = 2 photons 
and 9 = 0. 

FIG. 5. Cavity diffusion coefBcients versus atomic position, 
for (from the top curve going down) m = 8, 4, and 2 photons, 
with the atom-probe detuning A = 50 MHz (again G = 0). 

FIG. 6. (a) Simulated atomic trajectory and correspond- 
ing heterodyne signal for a single transit with m = 4, O = 0, 
and A = 50 MHz. Here x is atomic position along the cavity 
standing-wave (measured in units of the optical wavelength) 
and Vx is the atomic velocity along the standing-wave (mea- 
sured in optical wavelengths per /is). Note that the simulation 
includes three-dimensional classical center-of-mass motion for 
the atom, although only x and are displayed above. The 
two subplots in part (b) display the expectation value of the 
intracavity phase quadrature amplitude {qp) (which is pro- 
portional to Xp) filtered down to an analog bandwidth ~ 300 
kHz (left), as well as (qp) plus an appropriate amount of Gaus- 
sian noise to simulate shot noise in our heterodyne detection 
(right). 



FIG. 7. Simulated atom-transit signals (see text), dis- 
played at an analog bandwidth of 200 kHz. The atom-probe 

detuning A and probe power m are indicated for each sub- 
plot, and the atom-cavity detuning O is zero in all cases. 
Upper traces represent the intracavity amplitude quadrature 
(qa), with Gaussian noise added to reproduce the overall sig- 
nal-to-noise ratio predicted by equation (10) (yielding a quan- 
tity proportional to Xa)- Lower traces (with zero mean) rep- 
resent the intracavity phase quadrature (qp) plus Gaussian 
noise (proportional to Xp). 

FIG. 8. Individual atom-transit signals, displayed at the 
full analog bandwidth of 300 kHz and sampling rate of 10 MHz 
(12 bit resolution). The atom-probe detuning A and probe 
power m (see text) are indicated for each subplot, and the 
cavity-probe detuning O is zero in all cases. Upper traces rep- 
resent the amplitude quadrature Xa, whereas the lower traces 
(with zero mean) represent the phase quadrature Xp. Note 
that we have displaced the Xa traces by -1-400 in order to pre- 
vent them from overlapping with the Xp traces, and that the 
photocurrents are displayed in arbitrary units. 

FIG. 9. Dependence of transit phasor shapes on detuning 
(the values of probe power m and detuning A are indicated 
above each subplot). Each subplot displays an overlay of two 
data traces (gray spots), quantum- mechanical theory from 
the Master Equation (solid curve ending in a o, computed 
using equation (1)), and semiclassical theory (computed using 
equation (7), solid curve ending in an "x"). 

FIG. 10. Transit phasors for fixed detuning A = 10 MHz, 
with variable probe strength (as indicated). Each subplot 
displays an overlay of three data traces (gray spots), quan- 
tum-mechanical theory from the Master Equation (solid curve 
ending in a o, computed using equation (1)), and semiclassi- 
cal theory (computed using equation (7), solid curve ending 
in an "x"). 
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